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ABSTRACT 

Giant planets are usually thought to form within a few tens of AU of their host stars, and hence it came 
as a surprise when we found what appeared to be a planetary mass (~ 0.008 M ) companion around the 5 
Myr-old solar mass star 1RXS J160929. 1-210524 in the Upper Scorpius association. At the time, we took the 
object's membership in Upper Scorpius — established from near-infrared, H- and K-bwd spectroscopy — and 
its proximity (2.2", or 330 AU) to the primary as strong evidence for companionship, but could not verify their 
common proper motion. Here, we present follow-up astrometric measurements that confirm that the companion 
is indeed co-moving with the primary star, which we interpret as evidence that it is a truly bound planetary mass 
companion. We also present new /-band spectroscopy and 3.0-3.8 /im photometry of the companion. Based 
on a comparison with model spectra, these new measurements are consistent with the previous estimate of 
the companion effective temperature of 1800 ±200 K. We present a new estimate of the companion mass 
based on evolution models and the calculated bolometric luminosity of the companion; we obtain a value of 
0. 008^5'™! Mq, again consistent with our previous result. Finally, we present angular differential imaging 
observations of the system allowing us to rule out additional planets in the system more massive than 1 Mj up , 
2 M Jup and 8 M Jup at projected separations larger than 3" (-440 AU), 0.7" (-100 AU) and 0.35" (-50 AU), 
respectively. This companion is the least massive known to date at such a large orbital distance; it shows 
that objects in the planetary mass range exist at orbital separations of several hundred AU, posing a serious 
challenge for current formation models. 

Subject headings: stars: pre-main sequence — stars: low-mass, brown dwarfs — planetary systems 



1. INTRODUCTION 

The detection and characterization of low-mass substellar 
companions on large orbits (over several tens of AU) around 
stars is of great importance for our understanding of planet, 
brown dwarf and star formation, as well as for our understand- 
ing of the dynamical evolution of such companions in multi- 
ple systems or in circumstellar disks. Currently most forma- 
tion models of planets or low-mass brown dwarf companions 
- core accretion (e.g. Pollack et al. 1996), gravitational in- 
stability (e.g. Boss 1997) and fragmentation of a pre-stellar 
core (e.g. Padoan & Nordlund 2002; Bate et al. 2003) - can- 
not easily explain the existence of such companions at large 
orbital separations. A good determination of the semi-major 
axis, mass and mass ratio distributions of far out giant planets 
will help us understand these mechanisms in more detail and 
disentangle the role played by each one. 

Over the past decade, observations of young stars have 
led to the discovery of several companions with separations 
of several tens to hundreds of AU and masses only slightly 
above the —13 Mj up deuterium burning limit: TWA 5 B 
(Lowrance et al. 1999), AB Pic B (Chauvin et al. 2005b), 
GQ Lub B (Neuhauser et al. 2005), DH Tau B (Itoh et al. 
2005), CHXR 73 B (Luhman et al. 2006), LP 261-75 B (Reid 
& Walkowicz 2006), HN Peg B (Luhman et al. 2007), and 
CT Cha B (Schmidt et al. 2008). For higher mass brown 
dwarf companions, the list of known companions with sim- 
ilar orbital separations is even longer, the reader is referred to 
Zuckerman & Song (2009) for a recent compilation. More re- 
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cently, the mass limit of known substellar companions on very 
wide orbits has been pushed down below the deuterium burn- 
ing limit — i.e. in the planetary range — by a rapid succes- 
sion of discoveries, an —8 Mj up planetary mass object around 
the young solar analog 1RXS J160929.1-210524 (Lafreniere 
et al. 2008; hereafter called Paper I), three —7-10 Mj llp plan- 
ets around the A-type star HR 8799 (Marois et al. 2008), a 
— 1 Mj up planet orbiting inside the dust belt of the A-type star 
Fomalhaut (Kalas et al. 2008), and an —8 Mj up planet candi- 
date within the circumstellar disk of the A-type star f3 Pictoris 
(Lagrange et al. 2009). The newly-found companion to the 
Sun-like star GJ 758 (Thalmann et al. 2009), a 10^40 M Jup ob- 
ject orbiting at a separation of 29 AU, could add to the list of 
planetary mass companions if its true age turns out to be at the 
young end of current estimates. A few planetary mass objects 
have also been found in large orbits around brown dwarf pri- 
maries: 2MASS 1207b (Chauvin et al. 2005a), 2MASS 0441b 
(Todorov et al. 2010) and possibly UScoCTIO 108 b (Bejar 
et al. 2008). Beside actual planet detections, observations of 
asymmetries or ring-like structures in some debris disks offer 
indirect evidence for the presence of planets on large orbits 
around stars. Good examples include f3 Pic, HR 4796, e Eri, 
Vega, HD 141569, rj Corvi, Fomalhaut, and AB Aur (Hol- 
land et al. 1998; Greaves et al. 1998; Weinberger et al. 1999; 
Schneider et al. 1999; Kalas et al. 2005; Wyatt et al. 2005; 
Oppenheimer et al. 2008; see also Wyatt 2008). The pecu- 
liar disk morphologies observed in those systems could arise 
from gravitational perturbations by unseen planets. Modeling 
indicates that Neptune to Jupiter mass planets at several tens 
of AUs are typically needed to reproduce the observations. 
Although the current picture is still far from complete, all of 
these results offer a very interesting first glimpse of the full 
extent of orbital separations and masses of giant planets and 
low-mass brown dwarf companions around stars. 
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The planetary mass candidate companion to the young solar 
mass star 1RXS J 160929. 1-2 10524, which is a member of the 
5 Myr-old Upper Scorpius association (Preibisch et al. 2002) 
is the subject of this paper. The candidate companion, lying 
only 2.2"(^330 AU projected) from the primary star, appears 
to be a young planetary mass (0.008 M Q ) object in the Up- 
per Sco association, based on near-infrared imaging (JHK) 
and spectroscopy (HK). In particular, the triangular shape of 
its //-band spectrum clearly indicates low surface gravity, and 
the overall spectrum is in excellent agreement with the spectra 
of known low-mass brown dwarfs in the same association as 
well as with model spectra with T e ff ~ 1800 K and \ogg ~ 4, 
consistent with the inferred age and mass. While at the time of 
discovery, we had no proper motion measurements and hence 
could not confirm it was gravitationally bound, the proxim- 
ity and youth of the candidate companion offered very com- 
pelling evidence that the two form a bound pair. In this paper, 
we present follow-up astrometric observations that confirm 
the common proper motion of the candidate companion with 
the primary star. In addition, to improve the characterization 
of this system, we present new photometric observations at 
3.05 /im and 3.78 /m (L 1 ), new spectroscopic observations in 
the /-band, and new high-contrast angular differential imag- 
ing (ADI, Marois et al. 2006) observations. 

2. OBSERVATIONS AND DATA REDUCTION 

All observations used for this paper 3 were obtained with 
the Gemini North telescope and the ALTAIR adaptive optics 
system (Herriot et al. 2000). The target star itself was used for 
wave front sensing and the field lens of ALTAIR was used to 
reduce off-axis Strehl degradation due to anisoplanatism. 

2.1. Imaging 

All imaging observations were obtained with the NIRI cam- 
era (Hodapp et al. 2003) in //32 mode. Used with the AO 
system and the field lens, this yields a pixel scale of 21.4 mas 
and a field of view of 21 .9" x 21 .9" for the full detector array. 

The initial JHK imaging observations of 1 RXS J 1 609-2 1 05 
were obtained in 2008 April and June and are detailed in Pa- 
per I. For astrometric follow-up of the candidate companion, 
the star was observed again on 2009 April 6 and 2009 July 
1. These observations were made in the K' filter and used 
a 5-point dither pattern corresponding to the center and four 
corners of a square of side 10". At each dither position, we 
obtained one co-addition of thirty 0.3 s integrations in fast, 
high read-noise mode, and one 10 s integration in slow, low 
read-noise mode. At each position this provides an unsatu- 
rated image of the target star and a much deeper image of the 
field that can be spatially registered and scaled in flux without 
ambiguity. 

These imaging data were reduced using custom IDL rou- 
tines. A sky frame was constructed by taking the median of 
the images at all dither positions after masking out the re- 
gions dominated by the target's signal. After subtraction of 
this sky frame, the images were divided by a normalized flat- 
field. Then isolated bad pixels were replaced by the interpo- 
lated value of a third-order polynomial surface fit to the good 
pixels in a 7 x 7 pixels box while clustered bad pixels were 
simply masked out. Next the images were distortion corrected 

3 In addition to the observations described in this paper, we also obtained 
exploratory VLA observations at 21 and 6 cm as the NVSS catalog lists a 
radio source with a brightness of 5 mJy near the source. However the radio 
source is not associated; at 6 cm, the 0.7 mJy point source is at aj20O0 = 16 : 
09 : 29.6772 ± 0.001 1 and 5 J200 o = -21 : 04 : 49.591 ± 0.027. 



using the distortion solution provided by the Gemini staff. 4 
Finally, the long-exposure images were properly scaled in in- 
tensity and merged with their corresponding short-exposure 
images. 

Follow-up observations of the system were also obtained at 
3.05 ^m (H2O ice filter, 5% bandpass) and 3.78 ^m (Lf filter) 
on 2009 May 2 and 2009 June 8, respectively, to constrain 
the spectral energy distribution of the candidate companion 
further. For 3.05 /xm, 20 co-additions of ten 5 s integrations 
were obtained over 20 dither positions covering a square of 
side 8". To enable a shorter integration time and avoid satu- 
ration for 3.78 /im, only the central 768 x 768 pixels of the 
array was read and 30 co-additions of three hundred 0.2 s in- 
tegrations were obtained over 30 dither positions covering a 
square of side 6". 

For the long wavelength imaging, the data reduction in- 
volved first removing the sky background by subtracting, for 
each image, the mean of the two images acquired closest in 
time for which the star is not located in the vicinity of the 
candidate companion in the image considered. The images 
were then divided by a normalized flat field image created by 
taking the median of all images after masking out the star in 
each one and normalizing them to a median value of one. A 
~ 121 x 121 pixels median-filtered image was subsequently 
subtracted from each image to remove any low spatial fre- 
quency residual background signal, which was occasionally 
present; to avoid any potential bias in doing this, the median 
filter was computed after masking out a 9 x 9 pixel box over 
the star and companion. All images were finally aligned and 
combined using a 3a clipped average. 

An angular differential imaging (ADI, Marois et al. 2006) 
sequence of observations was obtained on 2009 April 6 to 
search for additional, closer and/or fainter companions in the 
system. These observations were done in the K' filter and 
only the central 512 x 512 pixel region of the detector was 
read. The main ADI sequence consisted of 320 exposures of 
10 s in low read noise mode; after every 20 of these images, 
two short, unsaturated co-additions of ten 0.3 s exposures in 
medium read noise mode were obtained to monitor the star 
position, the Strehl ratio, and the sky transparency. The total 
duration of the sequence is 1 .9 hours and the total field of view 
rotation is 39°. Both before and after the main ADI sequence, 
a 5-position dither pattern of short unsaturated observations 
(10 x 0.3s) was obtained to provide a sky frame for the reduc- 
tion of the short exposures as well as additional image of the 
point-spread-function (PSF) core. 

The short, unsaturated exposures for the ADI observations 
were reduced as above for the Zf-band imaging. The long, 
saturated exposures were dark subtracted, divided by a nor- 
malized flat field image, corrected for distortion using the so- 
lution provided by the Gemini staff, and padded with zeros 
to increase their size to 750 x 750 pixels and thus avoid loos- 
ing parts of the field of view during subsequent shifts and ro- 
tations. Bad pixels were corrected or masked as explained 
above. To compensate for possible variations in AO correc- 
tion and sky transparency, each long, saturated image was in- 
tensity scaled to a common flux in an aperture of radius 4 pix- 
els by interpolating the values measured for the short, unsat- 
urated images obtained through the ADI sequence; these in- 
tensity corrections were ~1 1% on average. The center of the 

4 The distortion correction used is centro-symmetric and is given by r = 
r' + 1.32 X 10~ 5 r' 2 , where r and r' are respectively the distortion-corrected 
and distorted radial pixel distances from the array center. 
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stellar PSF in all saturated images was found by maximiz- 
ing their cross-correlation with the closest unsaturated image 
in an annulus where both images are unsaturated; the images 
were then then shifted to position the star at the center of the 
frame. As a verification of the star positions determined in this 
way, we calculated the centroid of the planetary mass com- 
panion for all centered, de-rotated images and found a scatter 
of only ^0.12 pixels per axis. An azimuthally symmetric in- 
tensity profile was then subtracted from each image to remove 
the smooth seeing halo. Next, the stellar PSF speckles were 
removed from each image by subtracting an optimized refer- 
ence PSF image obtained using the "locally optimized com- 
bination of images" (LOCI) algorithm detailed in Lafreniere 
et al. (2007b). The individual residual images were then ro- 
tated to align their FOV and their median was taken to obtain 
the final residual image. 

2.2. Spectroscopy 

The spectroscopic observations in K and H, obtained on 
2008 June 21 and 2008 August 21 and 24 with NIRI, have 
already been presented in Lafreniere et al. (2009). The new J- 
band spectroscopic observations were obtained with the inte- 
gral field spectrograph NIFS (McGregor et al. 2003) on 2009 
July 3. The J grating was used with the ZJ blocking fil- 
ter, yielding a spectral resolving power of ^6000. Given the 
lower Strehl ratio in the J band and the proximity and con- 
trast of the candidate companion relative to the primary star, 
the observations were obtained in ADI mode to enable a sub- 
traction of the primary star PSF from every spectral slice of 
the data cube. The primary star was positioned near a cor- 
ner of the 3" x 3" NIFS field of view and the position angle 
of the instrument was adjusted to ensure that the candidate 
companion would remain within the field of view during the 
ADI sequence. The individual exposure time was 240 s in 
low read noise mode. The observations were split into two 
sub-sequences, of 15 and 10 exposures respectively, each one 
immediately followed by observation of an Ar lamp for wave- 
length calibration and of the AO star HD 155379 for telluric 
and instrumental transmission correction. Before the first se- 
quence and after the second sequence, two images were ob- 
tained with a dither of ±15" to provide a sky frame. A stan- 
dard NIFS flat field calibration sequence was also obtained. 

The reduction of the NIFS data, up to the reconstruction of 
the data cube, was made using the Gemini IRAF pipeline. The 
steps covered in this pipeline are sky background subtraction, 
flat-field and bad pixel correction, spatial and spectral cali- 
bration, and data cube reconstruction. The spatial sampling 
of the reconstructed data cube is 0.043" pixel" 1 . The instru- 
mental/telluric transmission correction and spectrum extrac- 
tion, which could be done using the Gemini pipeline, were 
instead done using custom IDL routines as described below. 
First the center of the stellar PSF were registered to a common 
position in all slices and all cubes of the sequence; the center 
positions were calculated by fitting a 2D Gaussian function. 
The stellar PSF was then subtracted, slice by slice, using a 
LOCI algorithm. The 'allowed' reference images included all 
those in which the companion had moved by at least 9 pixels 
(0.4") due to field of view rotation. The optimization region 
was a section of annulus with inner and outer radii of 1 .63" 
and 2.24" (14 pixel wide), respectively, excluding disks of di- 
ameter 8 pixels at each positions where the companion was 
present in any of the images to avoid any bias during the sub- 
traction. Figure 1 shows the collapsed cube before and after 
PSF subtraction. Following PSF subtraction, the spectrum of 
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Fig. 1. — Example of a median-collapsed NIFS data cube before (left) 
and after (right) PSF subtraction using ADI/LOCI as described in the text 
(§2.2). The field of view shown is 1.5" X 3.0" for each panel. The negative 
signals above and below the companion on the right panel are artifacts from 
the subtraction. 



the companion was extracted by summing the flux in a circu- 
lar aperture of diameter 5 pixels. The spectrum of the telluric 
standard was extracted using the same aperture, corrected for 
its spectral slope using a 9520 K blackbody curve, and its 
Paschen line at 1 .28 /xm was removed by dividing out a Voigt 
profile fit. The companion spectrum was divided by this spec- 
trum to correct for telluric and instrumental transmission. Fi- 
nally, the median of the 25 companion spectra was obtained. 

3. ANALYSIS AND RESULTS 

3.1. Verification of common proper motion 

The relative position and flux ratio of the star and compan- 
ion were determined in two ways which yielded values in very 
good agreement. First, the centroids were determined by fit- 
ting a 2D Gaussian function to both components and their flux 
ratio was found using aperture photometry. Second, the rel- 
ative position and relative flux were determined simultane- 
ously by minimizing the residual noise after subtraction of 
one component from the other. For a given observation, these 
measurements were repeated for all dither positions and their 
average was taken. The NIRI+ALTAIR with field lens pixel 
scale of 21.4 mas pixel" 1 and sky position angle from the im- 
ages FITS headers were used to convert image coordinates 
to separations and position angles. We did not observe as- 
trometric standards to monitor and calibrate possible changes 
in the pixel scale of the instrument between the epochs; such 
changes could occur, in principle, as a result of small changes 
in the shape or position of some components of the telescope, 
AO system or camera. Nevertheless, we can calculate up- 
per limits on pixel scale changes using several background 
sources in our field of view (see §3.3). For this purpose we 
used four background stars that have been well detected at 
all epochs. At each epoch, we computed the mean length 
and position angle of all baselines between these four stars (6 
baselines, lengths of 180-475 pixels). We then calculated the 
changes in values between the epochs for each baseline, and 
we finally computed, for each epoch, the average change over 
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TABLE 1 
ASTROMETRIC MEASUREMENTS 



Epoch 


Band 


P(") 


PA. (deg) 


2008.3205 


-Ks 


2.215 ±0.006 


27.75 ±0.10 


2008.4596 


H 


2.222 ±0.006 


27.76±0.10 


2008.4596 


J 


2.219 ±0.006 


27.76±0.10 


2009.2621 


K' 


2.222 ±0.006 


27.65 ±0.10 


2009.4972 


K' 


2.219 ±0.006 


27.74 ±0.10 



NOTE. — The position reported in the discovery pa- 
per was the mean of the first three rows. 



all baselines. The maximum pixel scale change observed was 
0.14%, corresponding to 3 mas at the separation of the com- 
panion (2.219" or 103.7 pixels), and the maximum change 
in position angle was 0.05°. These changes constitute up- 
per limits as they can be affected by non-zero proper motion 
of the background sources. Beside these possible changes in 
pixel scale, our measurements are affected by larger system- 
atic uncertainties which are estimated below. We also note 
that, as mentioned on the Altair webpage 5 , the plate scale 
of NIRI when used with the field lens of Altair is not accu- 
rately calibrated. This does not affect our analysis below but 
we mention it as it should be considered when comparing our 
measurements with data from other telescopes or other instru- 
ments. 

Table 1 provides a list of our astrometric measurements 
and figure 2 shows the changes in offset between the com- 
panion and primary over all epochs compared with the ex- 
pected changes for a distant background star, based on the 
proper motion and estimated distance of the primary; the fig- 
ure also shows the same measurements for four background 
stars. For the proper motion of the primary, we used the 
value given in the UCAC3 catalog (Zacharias et al. 2009): 
-11.2±1.5mas yr" 1 in RA and -21.9 ± 1.5 mas yr _1 in DEC. 
A few other values, differing by less than 6 mas yr" 1 per axis, 
are also available in the literature but have larger uncertain- 
ties; using those other values would not affect our conclu- 
sions. For the parallax motion we assumed the mean distance 
and spread of the Upper Sco association given in Preibisch 
et al. (2002), namely 145 ±20pc. As visible on the figure, the 
background star measurements are in good agreement with 
the background source model, although there are large sys- 
tematic errors for the individual data points. These systematic 
errors are likely due in part to residual geometric distortion. 
The average absolute difference between the measurements 
and the model for the background stars is 4.5 mas in both 
RA and DEC, or 6 mas in separation and 0. 1° in position an- 
gle; we adopt these values as our overall astrometric errors. 
Figure 2 also clearly shows that, within our uncertainties, the 
companion's position relative to the primary did not change 
over time; conversely, the measurements for the companion 
are inconsistent with the background model. From the earli- 
est to the latest epoch, the measurements are ~ 6ct off from the 
background source hypothesis. The \ 2 values for the common 
proper motion or background source hypotheses are 2.3 and 
33, respectively. A simple linear fit of the astrometry mea- 
surements with time yields changes of 2.0 ±5.6 mas yr for 
the separation and -0.05 ±0.09 deg yr -1 for the position angle 
(or -0.9 ±4.0 mas yr -1 in RA and 2.7 ±5.2 mas yr -1 in DEC), 
providing limits on the differential motion of the companion. 
Our precision is insufficient, however, to constrain the orbital 

5 http://www.gemini.edu/sciops/instruments/altair/field-lens-option 




Change in ARA (mas) 

FIG. 2. — Measured changes in offset between the companion candidate 
and the primary star (black filled circles) over the different epochs compared 
with the expected changes in offset for a distant stationary background source 
(solid black curve), as calculated from the proper motion and distance of the 
primary star. The measured changes in offset between four background stars 
and the primary star are also shown with small blue circles. The big red filled 
circles show the mean changes of the four background stars at each epoch. 
The blue and red circles are are connected with the background model curve 
at their corresponding epoch. The numbers inside the circles indicate the 
epoch according to the order of Table 1 . The error symbol at the lower left 
shows the uncertainty for individual data points. For all sources, the mean 
offset change over all epochs has been subtracted. 

motion of the companion, expected to be ^2.2 mas yr" 1 as- 
suming that the planet semi-major axis is close to its projected 
separation. 

Our new measurements clearly confirm the common proper 
motion of the two objects, but as pointed out in Paper I this 
does not in itself prove physical association as the velocity 
dispersion in the Upper Sco association is small. We previ- 
ously presented an analysis indicating that, conservatively, the 
probability that a planetary mass object in Upper Sco would 
fall within 2.5" of any of the stellar members we observed is 
only 0.002. We can now factor the common proper motion 
constraint into the analysis. Using the proper motion disper- 
sion of ^8 mas yr" 1 in RA and DEC reported by (Bouy & 
Martin 2009) for the low-mass members of Upper Sco, the 
probability that two unrelated Upper Sco objects would share 
a common proper motion within 4-5 mas yr" 1 per axis, as we 
have observed, is 0. 13. Thus taken together, the proximity and 
common proper motion of the primary and companion indi- 
cate a probability of 2.6 x 10" 4 of chance alignment. So given 
this very small probability, despite the lack of orbital motion 
detection, we treat the candidate companion as a truly bound 
companion for the remainder of this paper. 

3.2. Properties of the companion 

The new /-band spectrum of the companion is shown in 
Fig. 3 along with the spectrum of a field L3 dwarf, the spec- 
trum of a young L0 brown dwarf member of Upper Sco, and 
model spectra from the DUSTY (Chabrier et al. 2000) and 
DRIFT PHOENIX (Witte et al. 2009; Helling et al. 2008) 
atmosphere models for low and high surface gravity. The 
new /-band spectrum shows typical features of late-M or 
early-L dwarfs, namely important absorption by H2O beyond 
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FIG. 3.— 7-band spectrum of 1RXS J1609-2105 b (black) along with var- 
ious comparison spectra (red and blue). From the top, the first spectrum is 
compared with synthetic spectra from the DUSTY models with 7_j = 1 800 K 
and logg = 3.5 (red) or logg = 6.0 (blue); the second is compared with syn- 
thetic spectra from the DRIFT PHOENIX models with T cff = 1800 K and 
logg = 3.5 (red) or logg = 6.0 (blue); and the third is compared with the spec- 
trum of the LO brown dwarf USco J 160606-2335 13 (red), which is a known 
member of the Upper Sco association (Lodieu et al. 2008), and the spectrum 
of the field L3 dwarf 2MASS J1506544+132106 (blue, Cushing et al. 2005). 
The spectra are normalized at ~1.29 (im and binned to a resolving power of 
2000. 

~1.33 /im, the K I doublet at 1.24-1.25 /im, and absorption 
by FeH at 1 .24 /tm; the S/N is too small to identify other indi- 
vidual features but the general shape of the continuum is qual- 
itatively as expected through the band. As seen from the rea- 
sonable agreement with the model spectrum of T e s = 1800 K, 
the new spectrum is consistent with the previous estimate of 
the effective temperature of the companion. 

Figure 4 shows the merged JHK spectrum of 1RXS J 1609- 
2105 b along with synthetic spectra from the DUSTY and 
DRIFT PHOENIX atmosphere models for different temper- 
atures and surface gravities. From this figure, it is clear that 
the companion has low surface gravity: the slope of the con- 
tinuum through the H and K bands are much better fit by 
low gravity spectra for both sets of models. The DRIFT 
PHOENIX models provide a better fit of the overall spec- 
tral shape. In particular, the red side of the H band, which 
is poorly fit by the low-gravity DUSTY models owing to un- 
certainties in the model opacities used, is well fit the by the 
DRIFT PHOENIX models. For the latter models, the steep- 
ness of the slope on either side of the H band even provides 
some constraint on the effective temperature, favoring 1700- 
1800 K. Although less apparent from the figure, the K band is 
also better reproduced by the DRIFT PHOENIX models, the 
DUSTY models spectra falling off slightly too rapidly at the 
blue side. The shape of the K band spectrum, according to 
the DRIFT PHOENIX models, also favors a temperature of 
1800 K. 

The companion was detected in both the NB3.05 and L' 
images at the level of ^5cr and ~3<7, respectively; see Fig. 5. 
The flux of the companion and primary were determined from 
photometry in an aperture of diameter 6 pixels. The uncer- 
tainty on the companion flux was estimated from the dis- 
persion of the noise in 21 such apertures spread azimuthally 
around a circle of radius 2.2" centered on the primary star. 
We obtained a contrast of 6.8 ± 0.2 mag in NB3.05 and 
6.1 ±0.3 mag in L 1 . 

Unfortunately, we cannot derive the companion's flux in a 
straightforward manner as we did not observe photometric 
standard stars to calibrate the photometry, and no measure- 
ments of the primary flux in the same filters exist in the lit- 
erature. However there exist calibrated flux measurements of 
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FIG. 4.— Merged JHK spectrum of 1RXS J1609-2105 b (black) compared 
with various synthetic spectra from the DUSTY (left) and DRIFT PHOENIX 
(right) atmosphere models. From top to bottom, the comparison spectra have 
r eff of 2000 K, 1 800 K and 1 700 K. The red curves are for log g = 4.0 and the 
blue curves are for logg = 6.0. The spectra are normalized at 1.67 /im and 
binned to a resolving power of 500. Uncertainties in the colors of the compan- 
ion translate into ~15% uncertainty in the relative scaling of the companion 
spectrum between the bands. 
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FIG. 5.— Altair/NIRI images of 1RXS J1609-2105 in the 3.05 p,m nar- 
row band filter (left) and in the V filter (right). The white circle in each 
panel marks the expected position of the companion. The images have been 
convolved by a Gaussian kernel to improve the S/N of the companion. The 
negative regions left and right of the primary star in each panel are artifacts 
from the sky subtraction. 

the primary in several infrared bands that closely bracket the 
wavelength range of our observations: 2MASS J, H and K s 
(Cutri 2003), Spitzer IRAC 4.5 /tm and 8.0 /im, and Spitzer 
IRS peak-up imaging at 16 /im (Carpenter et al. 2006). We 
used these measurements, together with model spectra, to 
estimate the flux of the primary in the NB3.05 and L' fil- 
ters. More specifically, we used synthetic spectra from the 
NextGen models (Hauschildt et al. 1999) to compute syn- 
thetic fluxes in all the bands listed above, we then adjusted 
those synthetic fluxes to the measured values, and finally com- 
puted the synthetic average fluxes of the primary star in the 
NB3.05 and L' bands. For the 2MASS bands, the synthetic 
average flux densities were computed using the relative spec- 
tral response curves given on the 2MASS project webpage 6 
and the flux zero points determined by Rieke et al. (2008); 
for the Spitzer IRAC and IRS bands, the spectral response 
curves and zero points were taken from the respective instru- 
ments webpages 7 . Various synthetic spectra were adjusted 
to the photometric measurements by minimizing the reduced 
X 2 - The spectrum yielding the best fit, T e g = 4000 K and 
logg = 4.0 was used in combination with the spectral response 
curves of the NIRI NB3.05 and L' filters, taken from the NIRI 

6 http://www.ipac.caltech.edu/2mass/releases/allsky/doc/sec6_4a.html 

7 see http://ssc.spitzer.caltech.edu/irac/calibrationfiles 
and http://ssc.spitzer.caltech.edu/irs/calibrationfiles/ 
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website , to compute the primary star flux density in those 
two bands. We obtained values of 1.71 x 10~ 14 W m" 2 /j,m~ l 
for NB3.05 and 3.36 x 1(T 14 W m" 2 /urrT 1 for L'. Then 
using the observed contrast for the companion we obtained 
fluxes of (6.4 ± 1.3) x 10" 17 W m" 2 /an" 1 for NB3.05 and 
(6.2 ± 1 .8) x 1(T 17 W m" 2 fair 1 for L' . 

To verify the accuracy of the above procedure in estimating 
the primary star flux, we repeated it by ignoring one band (K s 
or IRAC 4.5 /im) when adjusting the synthetic fluxes to the 
observations and then retrieving the synthetic flux in that band 
from the adjusted spectrum. For both K s and IRAC 4.5 /im, 
the flux retrieved was within 3% from the measured value. 
Thus the error on the primary star flux resulting from this pro- 
cedure is a small contribution to the error on the companion 
flux. 

In Lafreniere et al. (2008) we derived a mass estimate for 
the companion primarily from our estimate of its effective 
temperature, and commented that the value found was consis- 
tent with the observed K s brightness. Here we adopt a differ- 
ent approach based on the bolometric luminosity of the com- 
panion inferred from a fit of model spectra to all photometric 
measurements available. The synthetic spectra were adjusted 
to the photometry in the same way as done above for the pri- 
mary. After adjustment to the data points, the model spectra 
were integrated over all wavelengths to obtain the total irra- 
diance, which was then converted to bolometric luminosity 
using the primary star distance estimate of 145 ± 20 pc. We 
repeated this procedure using both the DUSTY and DRIFT 
PHOENIX model spectra for a range of T e g from 1600 K to 
2000 K and \ogg from 3.5 to 6.0. For both sets of models, 
the spectra yielding the best agreement with the data points 
have T e ff =1700-1800 K and logg =3.5^4.5; representative 
spectra are shown in Fig. 6. Despite significant differences 
in the model spectra over the range explored, the range of 
bolometric luminosities obtained (±0.08 dex) is rather small, 
as noted by Marois et al. (2008) who used the same proce- 
dure for the HR 8799 planets. Then accounting for the un- 
certainty in distance, we obtain a bolometric luminosity of 
log(L/L Q ) = -3.55 ± 0.2. Compared to the prediction of the 
DUSTY evolution models (see Fig. 7), this luminosity indi- 
cates a mass of 0.008^' {J^ M©, in perfect agreement with the 
value reported in Lafreniere et al. (2008). For this mass, the 
dusty evolution models predict a T e ff of 1800 K, fully consis- 
tent with the best spectral fit. Thus, notwithstanding uncer- 
tainties in evolution models, the mass of 1RXS J 1609-2 105 b 
is in the planetary regime. 

We also used the same procedure, with NextGen model 
spectra of T eB =3800-4000 K (K7-M0) and logg =3.5-4.5, to 
estimate the bolometric luminosity of the primary star. We 
obtained a value of log(L/L ) = -0.37 ±0.15, in good agree- 
ment with the prediction (-0.36) of the evolution models of 
Baraffe et al. (1998) with a m [ x = 1.9 for an age of 5 Myr and 
Teff of 4000 K. The ratio of the companion bolometric lumi- 
nosity to that of the companion, which is independent of our 
distance uncertainty, is thus log(L comp /L p ,-i m ) = -3.18 ±0.09. 

3.3. Search for other companions 

The residual ADI image is shown in Fig. 8. In addition to 
the companion, this image shows 5 other point sources that 
are clearly detected. These sources are all further away from 
the primary and fainter than the companion. All of these five 
sources were detected at least twice in our multiple epoch 

8 http ://www. gemini . edu/sciops/instruments/niri/ 
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FIG. 6. — Photometric measurements of the companion (black points with 
eiTor bars) compared with the DUSTY model spectrum (red curve) for T e ff = 
1800 K and logg = 4.0 and the DRIFT PHOENIX model spectrum (blue 
curve) for T e ff = 1700 K and logg = 4.0. The short horizontal lines mark the 
corresponding model synthetic fluxes in each band; the length of the lines 
marks the extent of each photometric band. 




FIG. 7. — Calculated bolometric luminosity of the companion (point) 
compared with the DUSTY model evolutionary tracks for different masses 
(curves). 



TABLE 2 
Background sources detected 



P("r 


PA. (deg) a 


AK' (mag) 


4.240 ±0.006 


220.07 ±0.10 


8.40 ±0.10 


4.251 ±0.006 


347.57 ±0.10 


10.65 ±0.10 


4.687 ±0.006 


268.51 ±0.10 


9.50±0.10 


5.878 ±0.006 


46.40 ±0.10 


10.9 ±0.2 


6.706 ±0.006 


140.84 ±0.10 


11.4±0.2 



"At epoch 2009.2621. 

normal imaging and the astrometric measurements were suf- 
ficient to determine that they are background stars. For ref- 
erence, the position and contrast of these sources are given 
in Table 2. Visual inspection of the image did not reveal any 
other new, fainter point source. 

Detection limits in difference of magnitude as a function 
angular separation were computed from the standard devia- 
tion of the pixel values in an annulus of width equal to one 
PSF FWHM; the photometry was properly calibrated using 
the procedure described in §3.3 of Lafreniere et al. (2007a). 
The result is shown in Fig. 9. The evolution models of Baraffe 
et al. (2003), Chabrier et al. (2000) and Baraffe et al. (1998) 
were used to translate these detection limits into planetary 
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FIG. 8. — Residual ADI image normalized by the noise radial profile. The 
display is from -5a to +5a. The azimuthal positive/negative trails around the 
visible sources are artifact from the PSF subtraction algorithm. 
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FIG. 9. — Detection limits (5a) achieved by the ADI observations. The 
expected contrast for various planetary masses are shown toward the right of 
the figure. The filled circle marks the companion, 1RXS J 1609-2 105 b. 

masses, assuming an age of 5 Myr. Our results allow us to rule 
out the presence of planets more massive than 1 M Jup , 2 M Jup 
and 8 M Jup at projected separations larger than 3" (—440 AU), 
0.7" (-100 AU) and 0.35" (-50 AU), respectively. These 
limits indicate that if the planetary companion migrated out- 
ward (or was ejected) to its current location through planet- 
planet scattering, then the perturbing body, likely more mas- 
sive, would have to be located within 50 AU. 

4. CONCLUDING REMARKS 

The planetary mass companion 1RXS J 1609-2 105 b is the 
least massive known to date with an orbital separation of a 
few hundreds of AU. As discussed in Paper I, the existence of 
such a low mass companion orbiting so far away from its star 



poses a great challenge for star and planet formation models. 
Indeed, all main modes of formation -core accretion, grav- 
itational instability or binary-like formation- and migration 
-disk interaction or planet-planet interaction- face some ob- 
stacles when trying to account for such a companion. A lot of 
work has been done over the past two years in an attempt to 
explain the handful of companions at moderate-to-large sepa- 
rations recently found by direct imaging, the general outcome 
being that, although it remains challenging, given the right 
initial conditions, the above mechanisms can indeed succeed 
in producing distant gas giants. For example, gravitational 
scattering between multiple planets is a viable explanation for 
the planet we have found at —330 AU (e.g. Scharf & Menou 
2009; Veras et al. 2009), but of course this would require that 
several massive planets have successfully formed at smaller 
orbital separations (<30-100 AU). If this is the case, then 
1RXS J 1609-2 105 b should have a highly eccentric orbit and 
would likely eventually be ejected from the system. Alterna- 
tively, some authors have found that gravitational instability 
may also be a viable explanation, provided that the disk of this 
system was sufficiently massive, extended and had a relatively 
low opacity (e.g. Mem & Bate 2010; Boley et al. 2010); al- 
though admittedly, the regime where gravitational instability 
is effective remains a matter of debate (e.g. Kratter et al. 2010; 
Stamatellos & Whitworth 2008; Rafikov 2007; Boss 2006). 

More follow-up observations of this system over the com- 
ing years should help establish which, if any, of these mech- 
anisms can best explain the origin of this planet. In particu- 
lar, multi-epoch high-resolution spectroscopy of the primary 
and additional high-contrast imaging with improved sensi- 
tivity below 0.35" would allow to search for even closer-in 
companions; this would guide migration and ejection scenar- 
ios involving dynamical interaction between multiple planets. 
Also, higher S/N and/or higher resolution spectroscopy of the 
companion would allow to constrain the differential radial ve- 
locity of the primary and companion. This would constrain 
the orbital velocity of the companion -and its true physical 
bound- or its possible recent ejection. These data could also 
be used to constrain the metallicity of the planet, which could 
provide some clues to its formation mechanism. Better spec- 
troscopy would also enable a much more detailed comparison 
with model spectra. Along the same line, photometry over 
an extended spectral range would be useful to verify model 
spectra. A parallax measurement for the primary star would 
remove a big part of the uncertainty on the companion lu- 
minosity, which would then lead to a better mass estimate 
and/or a more constraining verification of evolution models. 
Finally, continued astrometry monitoring should begin to con- 
strain the orbital motion of the companion within a few years, 
given the expected differential motion of —2.2 mas yr" 1 . 
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